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Abstract 



A complete classification of integrable conservative hydrodynamic chains is pre- 
■ sented. These hydrodynamic chains are written via special coordinates - moments, 

lO ■ such that right hand sides of these infinite component systems depend linearly on 

| a discrete independent variable k. All variable coefficients of these hydrodynamic 

chains can be expressed via modular forms with respect to moment ^4°, via hy- 
pergeometric functions with respect to moment A 1 ; they depend polynomially on 
Q\ | moment A 2 and linearly on all other higher moments A k . A dispersionless Lax rep- 

resentation is found. Corresponding collisionless Boltzmann (Vlasov like kinetic) 
equation is derived. A Riemann mapping is constructed. A generating function of 
^ , conservation laws and commuting flows is presented. 
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1 Introduction 

In past years (2003 up to now) significant results were obtained in the theory of integrable 
hydro dynamic chains (see [13], [14], [6], [16], [19]). A first integrable hydrodynamic chain 

A\ = A k x +l + kA k ~ l A° x , = 0, 1,2, ... (1) 

was derived by D. Benney (see [2]) in 1973. An integrability of the Benney hydrodynamic 
chain can be illustrated by an existence of a generating function of conservation laws (see 

[2], [11]) 

where a generating function of conservation law densities is given by 

Ho_Hi_H : 
A A 2 A 



ii ill 112 /q\ 

P = X ~ — ~ — ~ T2 ~ •••> ( 3 ) 



whose all conservation law densities are polynomial functions with respect to moments A k , 
i.e. H Q = A , Hi = A 1 , H 2 = A 2 + (A ) 2 , H 3 = A 3 + A A 1 , ... It means that hydrodynamic 
chain (1) also can be written in the conservative form (see, for instance, [17]) 



— l—m I 

m=0 / . 



d t H = d x Hi, d t H k = ( H k+ i — - y H m H k _i_ m ) , k — 1,2, ... (4) 

We are interested in a description of integrable hydrodynamic chains written in the form 
(cf. (1)) 

A\ = hA k x +l + f A k + A k+1 (s A° x + siAl) + A k (r A x + nA x ) (5) 

+k{A k+1 (w A x + wiAl + w 2 A 2 x ) + A k (v A x + viA l x + v 2 A 2 x ) + A k -\u A x + mAl + u 2 A 2 x )], 

where coefficients fi,Sj,r k depend on first two moments A and A 1 only, while all other 
coefficients w m , v n , u p depend just on first three moments A , A 1 and A 2 . 

Recently, two particular cases of hydrodynamic chains (5) were completely investi- 
gated. The Hamiltonian hydrodynamic chains (here = dHi/dA k , k = 0, 1) 

A k = (a + (3)Hi,iA k+1 + PHi, A k x + [a(k + 1) + 2/3]A fc+1 (H 1 , 1 ) x + (ak + 2/5)A fc (H 1>0 ), 

are associated with the Kupershmidt Poisson brackets (see [5] and [10]); while the Hamil- 
tonian hydrodynamic chains (here H 2)k = dH 2 /dA k , k = 0, 1, 2) 



A k = 2H 2 , 2 ^ +1 + H 2>1 A k x + (k + 2)A k+ \^ 2 ) x + (k + l)A fe (H 2jl ), + kA k -\H 



2,0)x 



are associated with the Kupershmidt-Manin Poisson bracket (see the second part in [6], 
[8] and [10]). These hydrodynamic chains are integrable if and only if all components 
of corresponding Haantjes tensors vanish. It means that the corresponding Hamiltonian 
densities Hi (A , A 1 ) and H 2 (A°, A 1 , A 2 ) cannot be arbitrary. A full list of admissible 
expressions is given in [5] and [6], respectively. 

In a general case, the coefficient j\ in (5) is reducible to the unity, the coefficient f 
can be eliminate by an appropriate change of moments A k ; while all other coefficients can 
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be simplified in the integrable case only. Following the approach based on an existence 
of first three conservation laws and vanishing of the Haantjes tensor (see the first part in 
[6]), one can extract the integrable case 

A k = A k x +1 - k[(A k+1 + u A k + u_ 1 A fe - 1 )[ln(A 2 + a)] x - A k (u ) x - A k ~\u- X ) x ], (6) 

where functions Uq, it_i, a satisfy to an overdetermined system in an involution (see (54)). 
The same result can be obtained by the method of hydrodynamic reductions established 
by J. Gibbons and S.P. Tsarev in [9] and developed by E.V. Ferapontov and K.R. Khus- 
nutdinova in [1]. In this paper, we utilize the concept of the so-called symmetric hy- 
drodynamic reductions (see [16], [18]). In this case, an existence of a Riemann mapping 
X(q, A , A 1 , A 2 , ...) connecting the Vlasov type kinetic equation (see [7], [12], [21]) with 
hydrodynamic chain (6) leads to an overdetermined system in involution (36), (37), (38), 
(39), (40), (41), whose general solution can be parameterized by hypergeometric func- 
tions. Then a generating function of conservation laws can be found in quadratures. 
Thus, an infinite series of conservation laws densities H k {A° , A 1 , A k ) allows to rewrite 
hydrodynamic chain (6) in the conservative form 1 (cf. (4)) 

d t H k = d x F k (H ,H u ...,H k+1 ), fc = 0,l,2,... (7) 

We prove that its first two conservation laws coincide with first two conservation laws 
found in [6]. In the general case, E.V. Ferapontov and D.G. Marshall found that F 2 {Hq, Hi, H 2 ) = 
\n.H 2 + G(H , Hi) and functions F (H , Hi), G(H , Hi) satisfy to another overdetermined 
system in involution (see [6]). Moreover, we prove that this system in involution (55) 
is equivalent to system in involution (36), (37), (38), (39), (40), (41). Thus, a complete 
classification of integrable conservative hydrodynamic chains (7) is given in this paper. 

This paper is organized in the following way. In Section 2, symmetric 2N component 
hydrodynamic reductions are extracted by virtue of Zakharov's moment decomposition 
(see [18], [21]). Such 2N component hydrodynamic type systems contain N component 
symmetric sub-systems, which are still hydrodynamic reductions. These N component 
hydrodynamic type systems imply to the Vlasov type kinetic equation. We show that 
an existance of the Riemann mapping connecting this Vlasov type kinetic equation with 
hydrodynamic chain (6) allows to select all integrable hydrodynamic chains. In Section 3, 
canonical coordinates in a moment space are introduced. Then a further investigation sim- 
plifies. In Section 4, a special "triangular" case is completely integrated. Three variable 
coefficients in (6) can be parameterized by solutions of the so-called Halphen-Darboux 
system (see [1]). In Section 5, a generating function of conservation laws is found. In 
Conclusion, a generalization of the approach presented in this paper is discussed. 

2 Zakharov's moment decomposition 

The moment decomposition approach developed in [18] (see also [16]) is based on a concept 
of an existence of symmetric hydrodynamic type systems 

a l t = d x F(<L;p)\ p=a i, i = 1, 2, N, 

X A problem of a description of integrable hydrodynamic chains (7) was formulated in [13]. A particular 
and important Egorov's case F (H , Hi) = H i was investigated in [14]. 
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which are nothing but hydrodynamic reductions of the hydrodynamic chains, where in 
all known cases before (see, for instance, [5], [16], [17], [18], [19], [15]) each corresponding 
moment A k depends on N functions of a single variable 2 , i.e. 



N 



A k = J2fmk(a m ), k= 1,2,3,... 

m=0 

Another moment decomposition (introduced by V.E. Zakharov, see [21]) 



m=0 



also is applicable in all these known cases (see [18]). Benney hydrodynamic chain (1) under 
this moment decomposition reduces to the 2N component hydrodynamic type system 



\a°) , bt = (a i b%, 



2 

which possesses a formal reduction to the N component case (cf. (2)) 



if all field variables b k vanish, and A becomes a function of all rest field variables a n only. 
Let us replace a 1 by q(x, t, A), where A is a parameter. It means a 1 = q(x, t, £ l ), where £ l 
are arbitrary constants. Then (2) 

qt = qq x + A° x 

by a semi-hodograph transformation q(x, t, A) — > A(x, t, q) reduces to the linear equation 

A* = qX x — X q A° x1 

which is known as the Vlasov kinetic equation (see [21]; or the collisionless Boltzmann 
equation, see [7]). Suppose X(x,t,q) is a function X(q,A°, A 1 , ...), where all moments 
A k (x, t) satisfy Benney hydrodynamic chain (1). Since we suppose all moments A k are 
independent, one can obtain an infinite series of equations 

d k X = q~ k d X, A; = 0,1,2,..., (9) 

where d k = d/dA k , and (d q = d/dq) 

( ~ m Am-l\~ 1 

\ m=0 y / 



2 let us emphasize that N is an arbitrary natural number. 
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A solution of (9) is given by 3 



00 Am 
m=0 y 

where B\(q) and -82(9) are not determined yet functions. However, a substitution (11) 
into (10) yields Bi(q) = 1 and B 2 (q) = q. Then (11) becomes nothing else but an inverse 
series to (3). Thus, we conclude that a symptom of an integrability of hydrodynamic 
chains is an existence of a Riemann mapping \(q, A , A 1 , ...) connecting with Vlasov type 
kinetic equation (see below). In this paper, we utilize this property for a classification of 
integrable hydrodynamic chains. 

This moment decomposition approach can be extended on a wide class of hydrody- 
namic chains (cf. (5)) 

K M / K K \ 

A\ = fnA k x +n + E E Ak+Hs ™ + k E Ak+n ™n m A™ (12) 

n=0 m=0 \n=0 n=-l / 

where and M are arbitrary natural numbers, all functions fi,Sjk,wi p depend on first 
M + 1 moments (if A" = 1, M = 2, s 0>2 = 0, s 12 = and / , /1, s ,o, s ,i, si )0 , si ; i depend 
just on two first moments A , A 1 , these hydrodynamic chains reduce to (5)). Indeed, (12) 
reduces to N separate expressions for each index i (let remind that N is arbitrary) 

K 



n=0 

Af / K A" 

+ Yl E( a< ) fc+n6<s ™ + fc E ( ffii ) Hn ^- ) a; 



m=0 \n=0 n=— 1 



due to a substitution (8) in moments equipped by the index k only. Moreover, the above 
N expressions (due to their linear explicit dependence on a discrete variable k) can be 
split on two parts 



in 

X J 



K M K 

n=0 m=0 n=0 

K M K 

4 = E /» • (°t< + E E K) n+1 «w4™. (13) 

n=0 m=0n=-l 

As in the previous case, this 2N component hydrodynamic type system possesses N 
component reduction (13), where all moments A m and all variable coefficients depend 
on field variables a 11 only. 



3 It is well known that a general solution of the above linear equation is parameterized by one arbitrary 
function of a single variable A(A). However, in this approach, an existence of any solution is essential. 
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Our main observation successfully utilized in this approach is that integrable hy- 
drodynamic chain (12) is associated with the auxiliary equation 

K M K 

<?* = E fnQ n q x + E E ^"U" ( 14 ) 

n=0 m=0 n=-l 

which obtains due to a formal replacement a 1 — »• q in (13). It means, that equation (14) 
is compatible with hydrodynamic chain (12), where function q must depend on moments 
A k (x, t) and the parameter A. The semi-hodograph transformation q(x,t,X) <-» X(x,t,q) 
reduces (14) to the linear equation 

K UK 

h ~ E /«9" A * + EE ? re+1 ^nm^A 9 = 0, (15) 
n=0 m=0n=— 1 

which we call the Vlasov type kinetic equation (cf. [12]). The function X(x, t, q) depends on 
x, t implicitly via an explicit dependence on moments A k (x, t). We shall call hydrodynamic 
chain (12) integrable if a Riemann mapping X(q, A , A 1 , ...) connecting Vlasov type kinetic 
equation (15) with (12) exists. 

Examples: Hamiltonian hydrodynamic chains associated with the Kupershmidt- 
Manin Poisson bracket (see [8]; h n = dh/dA n , h nm = d 2 h/dA n dA m ) 

M-l M /M-l M-l \ 

A t = E^ + l ) h n+l^ +n + E E^ + l)^ fc+n ^+l,m + fc E ^/in+l.m A " 
n.=0 m=0 \n=0 n=-l / 

are connected with the Vlasov type kinetic equation 

M-l M M-l 

X t - E (" + ^^n+l^A, + E E 9 n+ ^n+l, m ^A g = 0, 
n=0 m=0n=— 1 

where the Hamiltonian is given by H = J h(A°, A 1 , A M )dx and (see (12)) 

fn = (n + l)h n +i, K = M - 1, S nm = (n + l)/l n +l,m, W«m = ^n+l,m- 

Hamiltonian hydrodynamic chains associated with the Kupershmidt Poisson brackets (see 
[5] and [10]) 

M M / M M \ 

A* = E(«« + £K^ +n + E E^ + 2/3)A fc+n /w + aA; E 

n=0 m=0 \n=0 n=0 / 

are connected with the Vlasov type kinetic equation 

M MM 

X t - E(«« + (3)h n q n X x + a E E <f +1 /w^A g = 

n=0 m=0 ri=0 

where the Hamiltonian is given by H = J h(A°, A 1 , A M )dx and (see (12)) 

f n = (an + ff)h n , K = M, s nm = (an + 2p)h nm , w nm = ah nm , w^ 1>m = 0. 
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Without loss of generality and for simplicity let us consider a hydrodynamic chain 
written in the form 



1 m / 1 1 \ 

A\ = fnA k x +n + E E Ak+US ™ + k E ^^nm A?. (16) 

n=0 m=0 \n=0 n= — l / 

If this hydrodynamic chain is integrable, then also all its higher commuting flows belong 
to the general class determined by (12) with appropriate choices natural numbers K and 
M. 

Lemma: The coefficient f\ can be fixed to the unity by the invertible point transfor- 
mation A k = (fi) k A k , then the coefficient f can be eliminated by the invertible point 
transformation 4 

k n,\ 

\k—m Ain 



m=0 ^ ' 



where is a binomial coefficient. If 8m f 1 = and 8m f° = 0, then hydrodynamic chain 
(16) reduces to the canonical form 



M 



A\ = A k+1 + E Ak+ns ™ + k E Ak+nw nm A™; (17) 



m=0 \n=0 



if djwf 1 7^ or 8m f° ^ 0, then (16) reduces to (17), but M replaces by M + 1, corre- 
spondingly. 

Proof: Hydrodynamic chain (16) is associated with the reduced version of (14) 

M 1 

qt = (fiq + fo)q x + E E <? n+1 ^ m A™. 

m=0 n=— 1 

Thus, the transformation q = fiq + f reduces the above equation to a more simple case 
with fi = 1 and f = 0. Corresponding Zakharov's moment decomposition (8) transforms 
accordingly 

N 

A k = Y,Ui^ + h) k b m . (18) 

m=0 

This is nothing else but a linear combination of aforementioned transformations. This 
point transformation A = A , A 1 = f\A° + (f ) 2 ,A 2 = (/i) 2 A 2 + + (/o) 3 , 

...,A M = (fi) M A M + M/o(/i) M ~ 1 ^4 A/ " 1 + - + (fo) M+ \ - cannot be inverted to a sim- 
ilar form due to complexity of functions fo{A°, A 1 , A ) and f\(A°, A 1 , A M ). Just 
higher moments A M+k (A°, A 1 , ...,A M+k ) became linear expressions with respect to higher 
moments A M+ \A M+ \... 

On the other hand, 2N component hydrodynamic type system (13) 

Ml M l 

b\ = Uia i +fo%+fibw x +b i E(°T s ^' < = (M 4 +/oK+E E («t +1 «w^ 

m=0 n=0 m=0 n=-l 



4 Similar transformations preserving the Kupershmidt-Manin Poisson bracket were considered in 
but with constant coefficients fo and f\. 



s 



under the aforementioned transformation c % = fid 1 + fo reduces to 

M 1 M+l 1 

b\ = C% + bV x + fi 4 X) J2( C TSnmA™, 4 = dc x + E E (cT +1 ™nmA™, (19) 
m=0 n=0 m=0 n=—l 

where 

Sim = — f dm m /l, SQ m = So m 7rS lm + fod m In /i — d m fo, 

h h 

Wim = - d m In /i, % M +1 = ty_i,M+i = /ltWo - /o^Af/l, 

Jl 

%m = W70m + (1 - 5 m , )5m-l/l - <9m/o + 2/ <9 m In / a - 2 



/l 

M I 1 1 \ 

X X A*+" Snm + X pA p+ ^ nm <9 P In / l7 

p=0 \n=0 n=-l / 



lU_l m = (l-^ m ,o)(/l9 m -l/o-/ 9 m -l/l) + /o9 m / -(/o) 2 g m In /!+ — „ Wlm -/oWOm + /lW-lm 

Jl 

M / 1 1 \ M / 1 1 \ 

-Jo E E + E ^ P+n ^- 3> 111 A+E E + E P AV+nW nm 9 p f . 

p=0 \n=0 n=-l / p=0 \n=0 n=-l / 

Due to (18), (19) can be written in the final form 

M 1 M+l 1 

b\ = c% + vj x + v J2( c T5n m A™, c\ = ee x + £ E ( C T +1 ^^ 



m=0 n=0 m=0 n=— 1 



where coefficients and tyjy are expressed via new moments A n . This is nothing else but 
a hydrodynamic reduction of hydrodynamic chain (17) 



M+l / 1 

A\ = A k x +1 + E E Ak+nS ™ + k E Ak+n ™nm K 



m=0 \n=0 ra=— 1 



Thus, Lemma is proved. 

Example: The remarkable Kupershmidt hydrodynamic chain (see [10], [15]) 

A\ = A k x +l + (3A°A k + (k + i)A k A° x , k = 0, 1, ... 

reduces to canonical form (17) 

i t fc = A k+1 + [(1 - p)h + 7 - /?]i fe i° + PkA*- 1 (A 1 + 7 ~^~ 1 (i°) 2 ^ , fc = 0, 1, ... 

Thus, we can investigate an integrability of hydrodynamic chain (16) written in a more 
convenient form (17) instead (16). In such (14) reduces to 

M 1 

ft^ + EE (f^nmA™. (20) 
m=0 n=— 1 



A consistency of (20) with (17) leads to an infinite set of equations 

d M+k q = q- k d M q, k = 0,1, 2, ... (21) 
and M + 1 equations (m — 0, 1, 2, M) reduces by virtue of (21) to 

AI-l / 1 1 \ 1 

(1-S mfi )d m ^q+Y^ ^s nm A k+n + w nm A k+n J d k q+E m -d M q = ^ w nm q n+1 +qd m q, 

k=0 \n=Q n=-l / ra=-l 

where M + 1 infinite sums are determined by 

s - = E (e s -^ M+n+fc + ( M + fc ) E ^«^ M+n+fc ) 4- ( 22 ) 



fc=0 \n=0 n=-l 
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However, all these infinite sums can be reduced to a sole sum only (see below). A consis- 
tency of the Vlasov type kinetic equation (cf. (15), see also (20)) 



M 1 



m=0 n=— 1 

with (17) yields an infinite set of equations (which is equivalent to (21) due to the trans- 
formation dkq = —dkX/dqX) 

d M +kX = q~ k d M X, k = 0,1, 2, 

whose solution is given by 

X = B 1 (q, A , A 1 , A M ~ l )\L + B 2 (q, A , A 1 , A M ~% 

where 

p=0 H 

while other M + 1 equations 

M-l / 1 i \ i 

(l-(5 m , )<9 m _iA+^ y^s nm A k+n + kJ2 w nm A k+n ) d k X+E m -d M X+ ^ u7 nm g n+1 9 g A = gd m A, 

fc=0 \n=0 n=-l / n=-l 

reduce to a linear system 5 G m (q, A , A 1 , A M )Y> + Q m (q, A , A 1 , A M ) = due to (22) 
expresses via (23) 

E Sn m q M+n+1 " E ( ra + l ) w n m q M+n A £ - £ ^ m g M+n+2 <9 9 £ 

v n=0 n= — 1 / n=— 1 

+ E U ' nm E n p-n-M ^ P ~ E Snm E 



g p-n-M / v nm / / gp-n-M ' 

n=-l p=0 ^ n=0 p=0 ^ 



5 This linear system does not contain a part proportional to d q Y>, because its corresponding coefficient 
vanishes authomatically. 
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Since both coefficients G m and Q m must vanish independently, a full 2M + 2 component 
system can be split on the two M + 1 component sub-systems of linear equations 



M-l / 1 



qd m In B l + (<5 m , - l)«9 m _! \nB 1 -J2[Yl s nmA k+n + k J] m„ m A fc+ " In £ x 

fc=0 \n=0 n=-l / 

1 1 

- w nmQ n+1 d q In Bi = ^[s nm - (n + l)w„ m ]g" (24) 



n=-l 



n=0 

and 



Af-l 



n=-l 



g^^+^o-l)^-!^- s nm A fc+ " + fc ^ u> nm ^ +n ) c\£? 2 - ^ w nm q n+1 d q B 2 

k=0 \n=0 n=-l 

1 

+ y^Km - (n + l)w nm ]q n B 2 = A°(s ln 



Wi r 



S 



m,0- 



n=0 



All derivatives of functions B\ and B 2 can be expressed from this linear system 
with variable coefficients and w^. A consistency of these derivatives leads to an 
overdetermined system in partial derivatives on Sy and Wki with respect to moments 
A , A 1 , A 1 and q, while a dependence on the highest moment A M can be found by 
a straightforward differentiation of linear system (24) written in the matrix 6 form 



* 



V 



q + * 

* 
* 



* —wi,iq — wo,iq + * 

* -wi j2 g 2 - w 0:2 q + * 

q + * -w h M-iq 2 ~ w ,M-iq + * 

* -w 1M q 2 - w M q + * 



\ / dolnB 1 
a In B 1 
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\ 



Om-i^B 1 
\dghiB 1 J 



( w q + * 
Wiq + * 

w M -iq + * 

\ w M q + * 



where Wk = s\k — 2w±k and the mark "*" means elements independent on q. Indeed, such 
a differential consequence is given by 



/ 


* 


* 


.. * 


-w' xx q 2 - w' 0l q + * 




* 


* 


.. * 


-w' h2 q 2 - w' 02 q + * 




* 


* 


.. * 




V 


* 


* 


.. * 


-w' 1>M q 2 - w' 0>M q + * 



\ / d \nB l \ 
d 1 In B 1 

Om-i^B 1 
J \ d.lnB 1 J 



( w' q + * \ 



w[q + * 



w'm-iQ + * 
\ ™' M q + * 



where the mark '"" means a partial derivative with respect to the moment A M . 
Lemma: Any row of the above linear system 

( *, *, *, -w' l m q 2 -w' m q + *, w' m q + *) 

is proportional to the last row from the previous linear system 



-w 1>M q 2 - w 0>M q + *, w M q + * ) 



(25) 



(26) 



6 a determinant of this (M + 1) x (M + 1) matrix is a polynomial of degree M + 2 with respect to q, 
except some special cases, like u>i,m = 0, which should be considered separately. 
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Proof: Indeed, let us consider the linear system 



/ q + * 


* 




-Wi^q 2 - w Qjl q + * 


— w q — * 


\ 


/ doB 1 




* 


q + * • 


. * 


-Wi^q 2 - w , 2 q + * 


—Wiq — * 




d 1 B 1 




* 


* 




- * -wi,M-iq 2 - w 0>M -iq H 


- * -w M -iq - 


* 


Om-iB 1 




* 


* 


. * 


-w 1)M q 2 - w 0jM q + * 


-w M q - * 




d q B x 




V * 


* 




-w'i, m <l 2 ~ w' ,mQ. + * 




/ 


V* 1 


I 



determined by the [M + 2) x (M + 2) matrix incorporating all rows of the original linear 
system and any row from its differential consequence. A determinant of this matrix equals 
zero for nontrivial solutions B 1 . Thus, the last row (see (25)) must be a linear combination 
of all other rows. However, most of them (m — 0, 1, M — 1) contain an element q + *, 
which does not exist in first M entries of this last row. Thus, the last row cannot be 
expressed via these higher flows except the row with the number M (see (26)). It means, 
that all elements of these two rows must be proportional to each other. Lemma is proved. 
Thus, the full set of equations is given by (n = 0, 1, M — 2) 



Pm 
Pm 



Pm 



m 



_ = 

hi 
JW-iy 

-TO / 

r t M 



c M-ly 
-M J 
M-l 



~-Mi 



M 



-, m = 0,1, ...,M- 1, 



(27) 
(28) 



M 



where 



Pr, 



S0,m — ^0,m + ( s l,m — 2Wi jm )q, 

nw_ hm A n ~ l + (s ,to + nw 0im )A n + (s 1>m + nw^ m )A 



ri+1 



8 m = w_ hm + w Q)m q + w hm q 2 . 
All these equations can be subsequently integrated. Indeed, the first ratio in (27) 

Pm 
Pm 



>M 

5m 



is nothing else but a cubic polynomial with respect to q. Since q is arbitrary, all four 
coefficients must vanish independently. A general solution of corresponding four ordinary 
differential equations (with respect to A M only) is given by 



So,m = (r - M + 1) w ltM , Si,m = (n - M + 1) wi 



Mi 



W ,M = U W hM , W- ljM = U-xW^m, (29) 

where functions r , Uq, r±, w_i depend on first M moments A , A 1 , A M ~ l . An integration 
of the second ratio in (27) 



Sm 



,M-ly 
, M-l 
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leads to 

wi,m = — ; , (30) 

a + A M ri 

where the function a depends on first M moments A , A 1 . . . , A 1 . An integration of 
the first ratio in (28) 

nl tr/ 
r'm U m 

leads to 

w- 1>m = u_iiOi, m + 7_ ljm , w , m = u w 1>m + 7 0m , (31) 

si, m = (n - M + 1) w lt7n + p 1>m , s , m = (r - M + 1) u; ljm + p m , 

where functions 7_i m ,7o m )Pom;Pim depend on first M moments A , A 1 , A M_1 . An 
integration of the second ratio in (28) 



A 1 i M-l 

Om 1 + e M 

leads to 

u m - A M p x 

Wi, m = — : , (32) 

a + A M n 

where functions ui m depend on first M moments A , A 1 , A . It is easy to see, that 
all other ratios in (27) and (28) are fulfilled by virtue of (29), (30), (31), (32). 

In the next Section, a more deep analysis is presented for hydrodynamic chain (5) 
written in form (17). 



3 Canonical variables 

The function B\ depends on first M moments A , A 1 , A M_1 only, but coefficients of linear 
system (24) depend also on A M explicitly via (29), (30), (31), (32). Thus, each derivative 
of In B\ can be expressed as a ratio of two polynomials with respect to q. In a general 
case (if W\ t M ^ 0), the common denominator is a polynomial of a degree M + 2. All 
numerators are polynomials of the same degree, except a numerator of derivative ln£>i 
with respect to q. Its degree is M + 1. Let us introduce roots qk{A°, A 1 , A M ~ X ) of this 
polynomial as basic field variables for further computations. In such a case, 

M+2 

^InB^-V;-^, (33) 
where a m (A° , A 1 , ...,yl M_1 ) are not yet determined functions. It means, that 

M+2 

Bi = a l[(q-q m )- a ™, (34) 

m=l 

where ao(A°, A 1 , A M ~ l ) is not yet determined function. A substitution (34) back to 
the first derivative of ln£>! with respect to q allows to express few (not all) variable 
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coefficients (r ,u ,ri,U-i,a, J_ 1>m , J 0>m , p 0>m , p 1>m , uj m , see the previous Section) via new 
field variables q k . Moreover, a$ and all other a m must be constant parameters (this is 
a consequence of an absence of logarithmic terms in derivatives of In B\ with respect to 
moments A k ); n is constant due to the constraint 

M+2 

Y> m = --M-1, (35) 

m=l 

following from comparison of r.h.s. in (33) with a corresponding expression from linear 
system (24). Without loss of generality, one can fix «o on the unity The compatibility 
conditions d k {d q hxBi) = d q (di-hxB 1 ),dk(d n lnB 1 ) = <9 n (<9fc lnSi) imply to explicit rela- 
tionships between some coefficients as well as dependencies d k q n via q m and rest of initial 
coefficients. Finally, the compatibility conditions dk(d m q n ) = dk(d m q n ) should lead to a 
parametrization of all coefficients via q m and their derivatives with respect to moments 
A°,A\...,A M ~\ 

However, this is not precisely true. Hydrodynamic chain (16) possesses a large class of 
invertible transformations, allowing to significantly reduce a number of distinguish coeffi- 
cients. For instance, hydrodynamic chain (5) contains 15 coefficients, while its integrable 
version (6) contains just 3 coefficients. It means, that transformation (18) is necessary 
but not sufficient for a most appropriate choice of reduced number of coefficients for a 
satisfactory investigation. To avoid complexity of this problem, in this Section, we restrict 
our consideration on the case M = 2 associated with hydrodynamic chain (5). 



4 General solution in the "triangular" case 

In this Section, we restrict our consideration on a most important case determined by 
the choice r = 1 and r\ = 1 (see (5) and comments to (12), i.e. the restrictions s ,2 = 
0, Si,2 = 0), i.e. (see (35)) 

4 

^2 a m = -2. 

m=l 

Moreover, we essentially can simplify further computations fixing all p kn = 0, where k,n — 
0,1. Nevertheless, this is not a particular case. A complete description of conservative 
integrable hydrodynamic chains (7) is given by (6). In general case (17), an infinite set of 
conservation laws can be written in the form (cf. (7)) 

d t H k = d x F k (H ,H 1 ,...,H M ), fc=0,l,2,...,M-l, 
dtHM+k = d x F M+k (H , Hi, H M+k+ i), k = 0, 1, 2, ... 

Just hydrodynamic chain (5) possesses an infinite set of conservation laws given by (7). 
Such hydrodynamic chains we call "triangular" in comparison with all other hydrodynamic 
chains (17), whose conservation laws possess a deviation from this triangular case, i.e. first 
M conservation law fluxes depend simultaneously on first M conservation law densities 
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As it was mentioned in the previous Section, the compatibility conditions d k {d q In B{) = 
d q (dk In Bi), dk{d n In Bi) = d n (d k lnBi) lead to the system in involution 

d iqk = ql + U ° q s k + U ~\ d 1 S=J2(2a m +l)q m , (36) 

m=l 



a m + 1 
q m 



(37) 



where 



1 / 4 44 

= — ( it_i y^(o; m + i)g m - JJgfe 

\ m=l fc=l m=l 

~ (g| + u g k + (a »g - u^) apM-i „ , , 

— o Co-uo, (do) 

= - J^(«fc + "m + l)q m qk, (39) 

m<fc 

<9 m-i = ( J^gfe - M-i ^(«fc + « m + l)gmg/t - («-i) 2 j , (40) 

\fc=l m<fc / 

4 

M = ^ a mqm, S = A°U-i + A^Q ~ (7, (41) 
m=l 

and 6 variable coefficients are connected with 3 others by (here k — 0, 1 only) 

7o,fc = ^fc^o, 7-l.fc = ^fcW-i; ^fc = -<9fc<x. 

In this case, all coefficients (29), (30), (31), (32) significantly reduce, then hydrodynamic 
chain (5) transforms to a more compact form given by (6). 
Let us introduce the auxiliary functions 

<lk - <?4 , , q 

q k = — - — , fc = l,2,3. 
Then equations (36) reduce to the form 

diq k = q k ^q k - + l)?mj , k = 1, 2, 3, (42) 

d^nS = ^2(2a m + l)q m , u_i = S I diq 4 - q 4 ^ a m9™ + ll- ( 43 ) 



m=l \ m=l 



A substitution u_i from the above system into (37) leads to the simple equation of the 
second order 

d\q± + (a 4 + l)gi<M 3 S = 0. (44) 
Let us introduce an intermediate function z = diq± and four functions c k (A°),k = 



0,1,2,3. 
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Lemma: A general solution of system (42) is given by 

q k = -di\n(z-c k ), fc = l,2,3, (45) 

where 

3 -2 3 

diz = c T\(z - Cm y- +1 , s = ^—T\(z - Cm y^ +1 \ (46) 

m=l m=l 

Proof: A substitution above formulas into (42), (43), (44) yields identities. 
A substitution (45) and (46) in (38) and (39) determines 

-1 3 



Q!4 + 1 A A , 

m=l 



where cq(A°) can be found by quadratures 



3 



(In c )' - - > y (a m + o 4 + 1) c„ 



m=l 

while all other Cfc(A fc ) satisfy a new modification of well-known generalized Darboux- 
Halphen system (see detail in [1]) 

. ax , «! + a 4 + 1 , 
Ci = — —r[ci [c 2 + c 3 ) - c 2 c 3 \ — q, 

«4 + 1 Oi4 + 1 

4 = — fr[ c 2 ( c i + c s) - cic 3 ] - " 2 + ft4 + _ c| ; (47) 

«4 + 1 a4 + l 

/ o 3 a 3 + a 4 + 1 2 

C 3 = — — T[ C 3 (Cl + C 2 ) - CiC 2 ] — — Cg. 

«4 + 1 Oi4 + 1 

Remark: A sole function 

S (A°) = (48) 
ci - c 3 

satisfies the so-called Schwarzian equation (see [1]) 

s'" 3(s") 2 _ ( ol\ol 2 + a 3 a 4 (a 2 + a 4 ) («i + a 3 ) (o 2 + o 3 ) (cci + o 4 ) ^ (s n2 



s' 2(s') 2 V s(s-l) s 2 (s-l) 2 

Then a solution of system (47) is given by 

Cl = -- ( In / . g ) ,02 = -- (ln^— ^ -) ,c 3 = --(\n(s a ^(s-lY 

Under the simple linear transformation 

ci = (1 + oji + a 3 )uj 1 - (1 + o 3 + a 4 )uj 2 + (1 - Oi + « 4 )u; 3 , 
c 2 = (1 - o 2 + 0:4)0; 1 - (1 + a 3 + a A )uj 2 + (1 + a 2 + a 3 )u 3 , 
c 3 = (1 + cci + 0:3)0; 1 + (1 - o 3 + o 4 )o; 2 + (1 + a 2 + o 3 )o; 3 , 
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the above formulas reduce to the form derived in [1], i.e. 

LOi = U 2 U 3 - U 1 (u 2 + U 3 ) + UJ 2 , 
U 2 = U\U 3 — U 2 (u\ + U 3 ) + U 2 , 

u 3 = u x u 2 - u 3 (ui + u 2 ) + u 2 , 

where 

u 2 = 0{(ui - u 2 )(u 3 - ux) + Pl(u 2 - u 3 )(u x - u 2 ) + Pl(u 3 - ux)(u 2 - u 3 ), 

Ul = ~\ ( ln i^i) ' " 2 = ~\ ( ln ^ny) ' " 3 = ~\ ( ln 7 

a 1 = ^(P 1 +P 2 -P 3 -l),a 2 = ^(-P 1 +P 2 + P 3 -l),a 3 = ^(P 1 -P 2 + P 3 -l) 
and (48) reduces to 



u 2 — u 3 

The function q^A , z) can be found by the quadrature 



3 -1 3 



dq, = c^z l[(z- c m T a ™~ x dz + Y[( z - cj— 1 P 2 {z)dA\ 

m=l m=l 

where P 2 {z) is a polynomial in z of the second degree, i.e. 

3 3 3 a + 1 3 
P 2 (z) = z 2 J^(a m + a 4 + l)c m + zjjc n — JJc n . 

m=l n=l m=l 171 n=l 

Then the first moment A 1 (A°,z) is determined by another quadrature 

dA 1 = c- 1 J] (2 - cj-^" 1 ^ + ( -*L_ J] (z - c m )- a "- 1 G 2 (^) + g 4 j dA°, 
where G 2 (z) is a polynomial in z of the second degree, i.e. 

3 3 3 

G 2 (z) = -z 2 + z ^2(a rn + a 4 + 2)c m + JJc n ^ 

m=l n=l m=l 

Thus, all functions uo,u_i and a (see (41), the second formula in (43) and the second 
formula in (46)) are expressed via above implicit dependencies z(A°, A 1 ), q^A , A 1 ), i.e. 

a = A°«_i + A'uo - -^frc* - c m y^^\ 

a A + 1 A \ 

m=l 



Or. 



U 



(«4 + l)c ^-C m A = l 



u_i 



7 2 n (* - c ^ _(2Qfc+i) + i l 4 ^ e n (* - _afc + <&■ 

(«4 + 1) (a 4 + 1) C ^ 2 - C m 



17 



5 Conservative hydrodynamic chains 

We believe that integrable hydrodynamic chain (6) must possess an infinite set of conser- 
vation laws (7). It means, that an infinite series of invertible triangular transformations 
Hk(A°, A 1 , A k ) can be found. In this Section, we present a generating function of 
conservation laws and prove an equivalence of the system in involution derived by E.V. 
Ferapontov and D.G. Marshall for conservative hydrodynamic chains (7) with system in 
involution (36)-(40) derived in the previous Section. 

Theorem: Integrable hydrodynamic chain (6) possesses a generating function of con- 
servation laws 

d tP (q,A 1 A 1 ) = d x Q(q J A° J A 1 ). (49) 

Proof: Integrable hydrodynamic chain (6) is associated with the Vlasov type kinetic 
equation (see (17)) 

Qt = QQx ~ (q 2 + u q + [log(A 2 + a)] x + q{u ) x + {u- X ) x . 

First two equations of hydrodynamic chain (6) are given by 

A° t = Al, A] = A 2 X - [{A 2 + uoA 1 + u^A°)[\n(A 2 + a)] x - A\u ) x - A (u^) x }. 

Since r.h.s. of q t and A\ contain derivative A 2 , formally Q should depend on A 2 . Differ- 
entiation (49) with respect to x and t 



d qP ■ q t + d oP ■ A° t + d lP ■ A] = d q Q ■ q x + d Q ■ Al + d x Q ■ A\ + d 2 Q ■ A 
leads to 7 

d q Q = qdgp, 

m fa , a (q 2 + qu + U- 1 )d o-\ 

doQ = d u-x + qd u — d g p 

\ a + A 2 J 

(a a a a (A 2 + A u^ l + A l u )d a\ 

V a + A 2 J 

in a , fa i a (g 2 + gw + ^-i)<V\ 

OiQ = O p+ diu-i + qd^o — d g p 

V ° + A 2 J 

a a a a (A 2 + A Q u_ l + A l u Q )d l a\ 

A)diu_i + AxdiUo — dtp, 

o + A 2 J 

o - AqU-i - AiUq q 2 + qu + u-! 

d iQ = —a d iP — 1 9 qP- 

a + A 2 cr + A 2 



If d 2 Q = 0, then 



q 2 + qu + u^ 

dip = - d q p 50 

a — AqU^i — AiU 



7 as usual, we are looking for a general solution. It means, that all A\ are considered independently. 
Thus, corresponding coefficients must vanish separately. 
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and all other above expressions simplify to the form 

a /-> a a /-> f(q 2 + qu + U-i) - d S) 
o q Q = qd q p, d Q = I h OqU-x + qd u ) d q p, 

an a , f (q 2 + quo + u-!) (u - diS) 

OiQ = 0oP+ 5 h + qdxuo d q p. 



A compatibility conditions di(d Q) = d (diQ), d\(d q Q) = d q (d±Q), d q (d Q) = d (d q Q) 
lead to three equations containing four second order derivatives do q p, d qq p, di q p, dooP only 
Taking into account (50), the derivative d\ q p is proportional to d qq p, and all other three 
derivatives d 0q p, d qq p, d 00 p can be expressed. Moreover, a direct further computation leads 
to the correspondence 

d qP = ^- (51) 

Thus, (see (34)) 

4 

d q p= Y[(q-q m ) am . (52) 

m=l 

The generating function of conservation law densities can be found in two quadratures 
(see (50) and (51)) 

4 2 4- 4- 4 

dp = J] (? - Qm)- am dq - q + ^ + M - 1 Y[(q- qJ'^dA 1 + (d p)dA°, 

m=l m=l 

where dop also is determined by corresponding second derivatives 

d(d p) = (d 0qP )dq + {d 01 p)dA l + (d 00 p)dA°. 

Nevertheless, an infinite series of conservation law densities can be found directly from 
(52). 

In contrary with the above approach, all conservation laws can be found iteratively. 
The zeroth conservation law is given by the zeroth equation 

d t H = d x F (H 0l H 1 ), 

such that A = H and A 1 = F (H , Hi) (see integrable hydrodynamic chain (6)). Let us 
introduce an intermediate notation h = doH\. 

Lemma: Integrable hydrodynamic chain (6) possesses first conservation law 

d t H x = d x \kH 2 + G{H Q ,H 1 )\, 

such that the second conservation law density 

H2= A^> 

the first conservation law density Hi can be found by two quadratures 



dH x = hdA° + ^dA 1 , dh = (d ^ - dA° + (d ^\ dA\ 



(53) 
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and the function G(H , Hi) is determined by the quadrature 

d S-u-i , (, , d 1 S-u 



dG = ^ s l dA° +{h+ ^—^Lj dA \ 
Proof: can be obtained by a straightforward computation. 

The system in involution on third derivatives of functions F(Hq, Hi) and G(Hq, Hi) 
was derived (see (55) in the Appendix) in paper [6]. 

Theorem: These functions F(H , Hi) and G(H , Hi) can be found in quadratures 

dG = ^oS-M-i + u h - Sh 2 ^j dH + (di In S - u + hs) dH u 

dF = -hSdH + SdHi, 
where d = dn ,di = d^- An inverse transformation is given by 

f) TP _ ft TPf) TP 

u = ^ d F - diG, S = diF, u„ x = d iF- . 1,1 + d Q Fd x G - &iFd G. 

diF diF 

Proof: can be obtained by a straightforward computation. 

Remark: All higher commuting flows belong to (12) in a general case. Indeed, a first 
commuting flow to hydrodynamic chain (6) is given by 

(A k+1 + A k UQ + A k ^i) - 

+^ [A k (u 0)Hl d x + d x u 0)Hl ) + A^ 1 (u- liHl d x + d x u- 1)Hl )] #2) 

where the third conservation law density is determined by (here Wo,#i = d^Uo, it-i,^ = 
d Hl U-i,a Hl = d Hl cr) 

1 3 

H 3 = (A 3 - au + A l u_i + -a Hl ) H% + -u Hl 

A compatibility condition (\t) y = {\)t of corresponding Vlasov type kinetic equations 
(15) (i.e. K — 1 and K = 2, respectively) leads to some 2+1 dimensional quasilinear 
equation of the second order (a general classification was presented in [3]), which will be 
considered in a separate paper. 



6 Egorov's case 

A most important and interesting case is the Egorov hydrodynamic chain (see [14]) se- 
lected by the simple choice Hi = A 1 (see (53)). In such a case, S = 1, then all = —1/2 
and general hydrodynamic chain (6) reduces to the form 

A^A k+1 -k[(A k+1 +A k diF+A k ~ 1 d F)[ln(A 2 +A 1 d 1 F+A°d F-l)} x -A k 
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where the function F is given by 



F = ^J V (A°)dA° + \n6 1 (A\A°). 
Here 77 (A ) is a solution of the Chazy equation 

7] = 67] — 2,7)7] 

and the Jacobi theta-function 

00 

e 1 (A 1 , A ) = 2 Y^(-l) n e~ in+1/2)2A ° sin[(2n + l)^ 1 ] 



n=0 



is connected with the above solution of the Chazy equation via an involutive system (see 
[14]) 



d 1 e 1 = -fie u d e 1 = ^ 2 -i)e l , 



dm = /, doti = ^J^P-AriP-8rfl-^rf'-^fjLl, 



d x l = \Ul 3 -4T]l 2 -8ri'l--ri", d l = I 2 - 7]l - r]' - -]i\l 4/ 3 - Ar]l 2 - 8rfl - -rf' . 
V 3 2 V 3 

7 Conclusion 

The crucial observation made in [18] is that a substitution of Zakharov moment decompo- 
sition (8) is applicable for hydrodynamic chains, whose r.h.s. expressions depend linearly 
on a discrete variable k and contain a finite number of common variable coefficients (see 
(12))- 

In comparison with approaches established earlier (see [4], [6], [9]), the method pre- 
sented in this paper is not universal but most effective. A complete classification of 
conservative hydrodynamic chains is given by virtue of their re-presentation in a special 
form (19). All conservation law densities H m can be expressed explicitly via moments A k ; 
all fluxes of corresponding conservation laws can be expressed explicitly via H m , ^11 com- 
muting flows can be constructed explicitly (and their conservation laws); infinitely many 
hydrodynamic reductions can be extracted. Thus, infinitely many particular solutions 
to integrable hydrodynamic chains (19) can be presented (by the generalized hodograph 
method, see [20]). 

8 Appendix 

The system in involution for functions «o, «-i and S describing a family of integrable 
hydrodynamic chains (6) possesses a general solution parameterized by 9 arbitrary con- 
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stants: 



a a a i e , o^- 1 ~~ u -i d i u o 
Oi iU-i = OiU-i ■ o± hi + 2 



di iUq = diu ■ d\ In 5" + 2 



£?o,iW-i = 9iit_i ■ <% In S — 2 
<9o i^o = 9\Uq ■ do In S + 2 



5 

<9iu_i + <9 u - u 9iM 

S ' 
u-i {diu-i + d u ) - u d u-i 



S 

S ' 



a c a a , d ° S ~ 2u -^ a c , u ° d ° S 
Oo,iS = <hu-i - o u H diS H — , 

dfS uodtS 2d S 
di,iS = — — + — — , (54) 

2d . S (ug - 2u_i) <9 S 9 M 5 
<9o,o£ = — ^ 1 ^ h u-i— ^— + - <9 w ) Mo 

+ (d u - diu-i - ^r<9 S' - u ~^ i ° ^ diS, 

a (a a \a «-i9i«o - 2d u_i 
»o,o«o = (<?o% - Oiu-i) »i«o H ^ <hS 

u-i (dxu-i + <9 u ) - u <9 u_i 2 (o>iit_i + d u ) - u diu 
~ 5 5 

a o« a /a a \a , «-i (di^-i + 2<9 u ) - 2 Uo<9 u_i 
d , M-i = 2d u_i • din - + d uo) o x u-i H dtb 

2(9 m„i - u^dxUo) + u Q diu^ u^dxuo - u- X (9d«-i + (#t«-i + <9 w ) Wo) + Wo^o^-i 
H ^ c ^+^ 

The system in involution describing conservative hydrodynamic chains (7) was derived 
in [6]: 

2%G (-4<9 G • dxG + 4d G ■ d F - d 2 F ^j d 0A G 



<%,o,oG 



diF d x F 



+ — =— 2 — ) 9i iG ~ 

d 1 F J ' d^o 

( _ 2 (d'G - 29 F ■ ^G + d 2 F + do,iF ) \ „ 
+ 2d G + =— ^ «9 , G, 

f<9iG - BoFq) 2 d ltl F _ 

di X iG = -± ^—L + Ad^G ■ d±F Q - d , F ■ d l F 

diF 

-2dxG ■ d 0jl F + 2d F ■ d 0A F - d G ■ d 1A F + (2 (d x G - d F ^ ^i,iG 
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<9 0)lil G = 2d G ■ B ltl G + 2<9 ,oG ■ d 1 F - 2d G ■ d 0il F 



(55) 



+ 



diF G 



d x F Q 



d ,o A G = -<9 2 G%^ -j 



diA.iFo 



d 1A F 



dhF 



dtFo 

^0,1,1-^0 = do nF ■ 8\Fq + OqG 



d 1 G - d F ) d 1;1 F + 2d 1 F ■ d 0>1 F , 
doiF 



diF Q 



di,iF , 



d F - d l G + %^ ) 



^0,0,0-^0 — d G 



d\G - 2d F ■ diG + d 2 F Q - d hl G 



doflFo 



dofiFo ~doG ■ d\G — 8q \G — 8qG ■ doF . 

2 = I d 0A F + 



dtFo 



d 2 G-d 0>0 G)d ljl F i 
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